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We present an overview of current-induced effects in nanoscale conductors with emphasis on 
their description at the atomic level. In particular, we discuss steady-state current fluctuations, 
current-induced forces, inelastic scattering and local heating. All of these properties are calculated 
in terms of single-particle wavefunctions computed using a scattering approach within the static 
density-functional theory of many-electron systems. Examples of current-induced effects in atomic 
and molecular wires will be given and comparison with experimental results will be provided when 
available. 
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CURRENT THROUGH A NANOSCALE 
JUNCTION 



Transport of electrical charge across a nanoscale junc- 
tion is accompanied by many effects, such as fluctuations 
of the average current; transfer of energy between elec- 
trons and ions and consequent heating of the junction; 
and forces on ions due to current-induced variations of 
the electronic distribution In this work we will dis- 
cuss these effects separately, and we will focus on their 
description at the atomic level. It is, however, important 
to realize that there has to be a (yet unknown) relation 
between these different properties. Such relation would 
constitute an important contribution to our understand- 
ing of transport in nanoscale systems. 

The static scattering approach to electrical conduc- 
tion will be the underlying theme of this review. How- 
ever, we point out that novel time-dependent formula- 
tions of the transport problem may lead us to a bet- 
ter understanding of these effects [2, 0] In what fol- 
lows we picture a nanoscale junction as formed by two 
semi-infinite electrodes held a fixed distance apart, with 
a nanoscale object bridging the gap between them. The 
nanoscale object could be a single atom, a chain of 
atomSja molecule, or any system with nanoscale dimen- 
sions lilEISSlElSEllilElIilElElES we 

then consider the problem of DC current flow from one 
electrode to the other as the result of an applied bias Vb ■ 
If we take the left electrode to be positively biased, elec- 
trons will flow from the right electrode to the left one. A 
buildup of negative charge will be present on the surface 
of the right electrode within a screening length of the 
electrode surface. Similarly, a buildup of positive charge 
will exist on the surface of the left electrode due to a 
corresponding depletion of electrons 0, 0] . We will as- 
sume that, far away from the junction, the electrons in 
each electrode are in local thermal equilibrium and their 
statistics are described by the Fermi-Dirac distribution, 
so that Vb = Efr — -E-FL, where -Efl(r) is the chemical 
potential deep in the left (right) electrode. (Here, as in 
the rest of this work, we use atomic units.) 

We calculate the transport properties of this system 
by expanding the stationary states of the Hamiltonian 



into a set of left- and right-moving waves. We then sum 
each left- and right-moving state, weighting them with a 
Fermi function according to their energy. The stationary 
scattering states of the bare electrodes have the form 



* £ K ll (r)= e iK "- Y n B K ll W ! 
with the following boundary conditions [l7|: 
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K|| is the electron momentum in the plane parallel to the 
electrode surfaces. We have defined hk R — E — i |K|| | 2 — 
u e //(oo), and \k\^ = E - | 2 - w e //(-oo). Y is the 
component of the position vector in a plane parallel to 
the electrode surfaces, and z is the coordinate perpendic- 
ular to them. v e ff (±oo) is the bottom of the electronic 
energy band deep within the right/left electrode. The 
wavefunctions 4"° satisfy the continuum normalization 
condition 



dM*W(r)r 



■EKn 



{r)=5(E'-E)S(K [l / -K ll ). (3) 



Since the details of the electrodes are not important 
up to the interface with the sample, we represent them 
using a uniform-background (jellium) model |14lll5l |. The 
potential V the electrons exp_erience when they scatter 
through the nano junction is 



V(r,r') = v ps (r,r') + 



3 ,; Mr") 



d-V 



r — r 



u xc [n(r)] - u xc [n°(r)] 

6(r-r>). (4) 



The term 



is the electron-ion interaction potential 



that we represent with (nonlocal) pseudopotentials; v xc 
is the exchange-correlation potential computed using the 
local-density approximation to density functional theory 
(DFT) j63; n° (r) is the electronic density for the pair 
of biased bare electrodes; n(r) is the electronic density 
for the total system, and <5n(r) is their difference. We 
are implicitely assuming that static DFT gives a reason- 
able account of the scattering properties of a nanoscale 
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system, at least in linear response. While this may be 
true for metallic junctions it is not obvious for molecular 
junctions |l8|. However, in linear response and far from 
the resonant regime, we expect basic physical trends for 
these systems to be reproduced well by static DFT [lij. 

The full Hamiltonian of the system is H = Hq + V, 
where Ho is the Hamiltonian due to the bare biased elec- 
trodes, and V is the scattering potential. Our next task 
is to find the self-consistent solutions to the equation 
H^e — E^Ei which we can put into the Lippmann- 
Schwinger form: 

*£K,i ( r ) = *ek„ ( r ) 

+ J d 3 r'd 3 r"G° E (r,r')V(r',r")^ EKll {r"), (5) 

The quantity G E is the Green's function for the bare 
electrodes, and needs to be calculated for each energy E. 
Lastly, the total density of the system is 



n(r) 



l*<ttl 



2 dE /d 2 K|||* BK|l (r)| 2 , (6) 



where we have included a factor of 2 due to spin degen- 
eracy. The \I/j's are the bound states of H, if any exist. 
They can be calculated by direct diagonalization of the 
full Hamiltonian H . In order to find a self-consistent 
solution for the density, equations J3J), J5J and © are 
solved iteratively. 

One has the freedom to choose a basis set to represent 
the wave functions. Due to the non-variational nature 
of the electrical current, the issue of which basis set to 
use in a given calculation is far from trivial [19(. In the 
examples that follow, plane waves have been used as basis 
set |61|, which allows easy testing of the convergence of 
the results 

Once the wavefunctions have been calculated self- 
consistently, the total electric current density (at zero 
temperature) is given by 

j(r)=-2 f ™dE I "d 2 K||Im{[* EK|l (r)]*V* BK|l (r)}. 

(7) 

The quantity we are interested in is the extra current 5 J 
due to presence of the nanojunction [ijj, El H3, E3| (two 
semi-infinite contacts have infinite surface area, and thus 
pass infinite current) 



5 J 



d 2 Yi-[j(r)-j°(r)], 



(8) 



where j is the electric current density, and we have de- 
fined j° to be the current density in the absence of the 
nanojunction. The current defined in equation (JHJ is the 
average current that flows across the nanojunction. Fluc- 
tuations with respect to that average are expected and 
will be described later in this work. 



It is worth noting that, in contrast to references [24LI2E 
EH. , we do not calculate the electrical current using 
any of the Landauer formulas [Ii||2j|; rather, the cur- 
rent is calculated as the expectation value of the current 
operator over single-particle states. Clearly, it can be 
proven that when the transmission probabilities are ex- 
tracted from the single-particle scattering wavefunctions, 
we recover from equation JSJ the conventional two-probe 
Landauer formula |3f)j |. 



CURRENT-INDUCED FORCES 

Now that we have computed the stationary scatter- 
ing states and corresponding current we can look at one 
of the effects induced by electron flow. We first focus 
on current-induced forces, i.e. the phenomenon by which 
atoms in a current-carrying wire are subject to forces due 
to the local change in self-consistent electronic distribu- 
tion [UII3 

There are several open questions related to 
current-induced forces, the most notable of which is their 
conservative nature |33j. Here we will focus on their de- 
pendence on some microscopic properties. We will show 
that these forces are a non-linear function of the junction 
properties, such as its current density, charge density and 
length. 
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FIG. 1: Different contributions to the total force for each of 
the three atoms composing a silicon wire, as a function of 
bias. See text for details. Reprinted figure with permission 
from |34|. Copyright 2003 by the American Physical Society. 

Let us first define forces in a current-carrying wire. 
We first note that in this non-equilibrium problem the 
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usual Hellman-Feynmann theorem is not valid |3lL 13; 
A meaningful definition of force on ions can then be ob- 
tained by either the classical limit of Ehrenfest's theorem 
applied to the rate of change of ionic momentum |3l| or 
from the Euler-Lagrange equation of motion for the clas- 
sical ions [35| . Both approaches yield the same expression 
for the force on an ion, with position R, due to the self- 
consistent electronic density p(r) under current flow [sEl , 



j dv ps 



(9) 



which can be equivalently written as 



■E 



dv 



ps 
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lim 



dE{ V>a 



dv 



ps 



<9R 



A 



(10) 

The first term on the RHS of Eq. IjlOl) is similar to the 
usual Hellmann-Feynman contribution to the force due to 
localized electronic states \tpi). The second term is the 
contribution due to continuum states |3l| . The wavefunc- 
tions IV'a) are eigendifferentials for each energy interval 
A in the continuum a [31j ■ Finally, an additional contri- 
bution from ion-ion interactions needs to be taken into 
account. 

For convenience we can separate the total force into 
two contributions. We define "direct force" to be the 
portion of the total force due to the charge distribution 
of the biased bare electrodes. We then call the remain- 
ing term the "electron wind force" [H^. There exists an 
inelastic portion of the electron wind force due to energy 
transfer from electrons to ions, in the form of phonon 
excitations; however, this contribution is generally negli- 
gible compared to the elastic one and will be neglected 
here 0. 

We note that atomic relaxations induced by current 
flow do not have a large effect on the absolute value of 
the current that passes through the junction, even in the 
limit of high voltages and current densities 0, 0, |3|| 
|33 • While this result says nothing about the mechanical 
stability of current-carrying wires under bias, it allows 
us to study several transport properties assuming atomic 
positions fixed at their zero bias value. 

Let us now look at the dependence of these forces on 
microscopic details. As an example, in Figure \I]we con- 
sider a nanoscopic wire composed of three silicon atoms 
between two bulk electrodes. We keep atomic positions 
fixed at their zero bias values. For each atom in the junc- 
tion, we plot the direct force, the electron wind force due 
to the continuum states, and the electron wind force due 
to the discrete part of the spectrum. 

The direct force is almost linear with applied bias, with 
the force on the central atom being the largest. Closer to 
the electrode surfaces, the electrostatic potential ceases 
to be a linear function of position, causing the force on 
the other two atoms to be slightly smaller. Note that the 
force on the central atom due to states in the continuum 



is positive (i.e. pushes the atom to the left, in the same 
direction as the electron flow). However, the force on 
the central atom due to states in the discrete spectrum 
is almost zero, even at larger biases. The reason for this 
behavior is that the amount of charge localized around 
the central atom is almost constant in the applied bias, 
and so the corresponding force is also constant |34|. 

Figure Ufa) through (c) shows the force on specific 
atoms in silicon wires composed of varying numbers of 
atoms, while Figure Efd) shows the average force on the 
total wires. It is difficult to extract overall trends from 
Figures Ufa) through (c), although it is interesting to 
note that the second atom from the left (labeled Si2) is 
consistently the atom that experiences the greatest force. 
Figure Eld), however, follows an obvious trend: the av- 
erage force saturates with increasing number of atoms. 
This comes from the fact that, as the length of the wire 
increases, the boundary effects due to the electrodes be- 
come less important. This result suggests that longer 




FIG. 2: Total force on each atom in atomic wires consisting 
of (a)two, (b)three, and (c) four silicon atoms. The inset 
shows a schematic of the wire that is composed of two silicon 
atoms, (d) The avera ge fo rce on each wire. Reprinted figure 
with permission from [341 . Copyright 2003 by the American 
Physical Society. 
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wires are more difficult to break under current flow, an 
important factor to consider in nanoscale devices. 



SHOT NOISE 

In this section we will discuss steady-state current fluc- 
tuations that occur due to the quantization of charge. 
Shot noise is quite distinct from (equilibrium) thermal 
noise and is generally unavoidable, even at zero temper- 
ature 38]. 

We will first derive an expression of shot noise in terms 
of single-particle wavefunctions. In order to do so, let us 
write the field operator for electrons \& as combination of 
a field operator due to electrons incident from the left, 
and one due to electrons incident from the right pffll l39l ] 



(11) 



We can further expand each of these field operators into 
the single particle wave functions JSJl: 



L(R) _ 



-iwt„L(R)^L(R) 



"E 



E 



'(r,K||). 



(12) 



Again, we have used K|| to denote the component of the 
electron's incident momentum that is parallel to the sur- 
faces of the electrodes 63] . 

The coefficients are the annihilation operators 

for electrons incident from the left (right) reservoir. 
These operators satisfy the usual anticommutation re- 
lation {d£, a^t} = 5$ X 5(E - E'), where 0, x = R, L. We 
will again assume that the electrons coming from the 
left (right) reservoir are in local thermal equilibrium at 
a temperature T e far away from the junction, so that 
their statistics are given by the Fermi-Dirac distribution 
function / 



L(R) 



i.e. 



6^ X 6(E - E')j% 
S^ X S(E - E') 

e [E-E F4> ]/k B T a _|_ i ■ 



(13) 
(14) 



Using the field operator 1)11(1. we can define the current 
operator 



I(z,t) = -i j dY J dK\\(¥d z i> -d z &i>). 



(15) 



In the limit of zero temperature, the Fermi distribution 
reduces to a step function; if we again assume the chemi- 
cal potential in the right reservoir Eyh is higher than the 
chemical potential in the left reservoir -Efl, then at zero 
temperature the average value of the current is just 



(I) 



dE dY 



dK \\lE.E> 



where 



f<P<x _ 
1 e,e' — 



(16) 



(17) 



and * denotes complex conjugation. 

We define shot noise as the Fourier transform of the 
electric current autocorrelation function in the limit of 
zero frequency and zero temperature |38l Eflj : 



2ttS(lu)= dte iwt (AI(t)AI(0)), 



(18) 



where the excess current operator A7(t) is equal to 
I(t) — (I). The evaluation of this expression is quite 
cumbersome and therefore we do not include it here. We 
only mention that we need to evaluate terms of the form 
(A4A3A2A1), where the A^s are raising and lowering op- 
erators. We can evaluate these terms with the Bloch-De 
Dominicis theorem ! 1 if : 

(AiAsAsAi) = (A 4 A 3 )(A 2 A 1 ) + V (A 4 A 2 }{A 3 A 1 ) 

+r) 2 (Aai)(A 3 A 2 ), (19) 

where rj = — 1 for fermions, and +1 for bosons. The end 
result is 

S(u>) = A^/l+J 1 -/!) /dYi/dK! 

*I$U,,e Ay 2 [dKjx* (20) 



In the limit of zero frequency and zero temperature, 
equation l(2T)|) reduces to pij, |3j| 



S = 



dE\ 



-E-FL 



dR / dKI^ K E 



(21) 



This is the desired expression relating shot noise to single- 
particle wavefunctions [frij] . 

In the case of uncorrelated electrons, the magnitude of 
the shot noise becomes Sp = 2el, where e is the electron 
charge, and I is the dc current. This corresponds to the 
noise of a series of incident particles, the time between the 
arrivals of which follows a Poissonian distribution func- 
tion. For this reason, 5*p is sometimes called the Poisson 
value for the shot noise In general, however, the 

magnitude of the shot noise will be less than the Poisson 
value. A relative measure of noise is therefore the Fano 
factor F, defined as the ratio between shot noise S and 
the Poisson limit Sp. 

We are now ready to discuss an example of noise prop- 
erties of a nanoscale junction. We again look at the prop- 
erties of a junction formed by a short wire of silicon atoms 
between two bulk electrodes. Other examples of noise in 
atomic-scale systems can be found in references [3{J and 

Figure [3] shows the results of conductance and Fano 
factor for such a system for a bias of 0.01 V. We first 
notice that the Fano factor is strongly nonlinear as a 
function of bias. It is also considerably enhanced for very 
short wires due to the large contribution from the metal 
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FIG. 3: Fano factor and conductance for a nanojunction com- 
posed of different numbers of silicon atoms. Reprinted figure 
with permission from 39]. Copyright 2003 by the American 
Physical Society. 



electrodes. In addition, the Fano factor oscillates as a 
function of number of silicon atoms. The conductance 
shows a similar but opposite oscillatory trend. Both the 
noise and conductance oscillations are due to the fact 
that a Si wire made up of an even number of atoms has 
fully-occupied it orbitals, while a wire made up of an odd 
number of atoms has a half-filled it state at the Fermi 
level 0. 

Other interesting properties of shot noise in nanoscale 
systems include its dependence on contact geometry and 
interwire interactions [3(J ■ We refer the reader to the 
original papers for a discussion of these effects. 



LOCAL HEATING 

Local heating occurs when electrons in a current- 
carrying wire exchange energy with phonons [lil liH litj . 
I47I |4S| . Accordingly, there are four main processes that 
contribute to local heating in the junction: 

1. Cooling processes in which an electron incident 
from the left absorbs a phonon, which we will de- 
note by the superscript L, 1; 

2. Heating processes in which an electron incident 
from the left emits a phonon, denoted by the su- 
perscript L, 2; 

3. Cooling processes in which an electron incident 
from the right absorbs a phonon, denoted by the 
superscipt R, 1; and 

4. Heating processes in which an electron incident 
from the right emits a phonon, denoted by the su- 
perscript R, 2. 

The power generated in a nanoscale junction is given by 
the sum of the average power generated by those four 



{{W R ' 2 ) + (W L ' 2 ) - (H^ 1 ) - {W^ 1 )) 



vib.mo 



clcs 



(22) 

In addition to the above processes, cooling occurs due 
to dissipation of energy in the bulk electrodes [S^. In 
order to evaluate the power generated by each of those 
processes, we first need to consider the full many-body 
Hamiltonian of the system: 



H — H c \ + H ion + H c 



(23) 



Here, H e \ is the electronic Hamiltonian (including 
electron-electron effects), Hi on is the ionic Hamiltonian, 
given by 

N p2 

H ™ = J2ut + T,Kon(&i-B.j), (24) 



and H c ]^\ nn is the electron-ion interaction, 



ion — ^ ^ ^cl— ion(l"i Rj i) 

id 



(25) 



Pj , Mj and R^ denote the momentum, mass and position 
of the i th ion (out of a total of N ions), while r$ denotes 
the position of the i th electron. 

We assume that each ion executes a small vibration 
about its equilibrium position so that its diplacement 



is given by Qj 



R , 



R'J. We can decouple these ionic 



vibrations by introducing normal coordinates {qjp} so 
that the a th component (a = x, y, z) of Q,; is given by 



N 3 

i=i 0-1 



(26) 



The coefficients A ia jp obey the orthonormality rela- 
tions Y,i, a M i A iot,ol3 xA ia,j'p> = ^/3j'/3'- We now can 
solve this problem in the usual way by introducing bo- 
son creation and annihilation operators b-p and bjp for 
the (j/?) th mode; the creation and annihilation opera- 
tors obey the commutation relation [bjf3,b^ l/3 ,] ~ Sjp^'p/. 
With this transformation, the total ionic Hamiltonian is 
just the sum of the Hamiltonians for each normal mode: 



(27) 



Next, much like in Section D we expand the field oper- 
ator for the electrons into a part that describes electrons 
incident from the left, and a part that describes electrons 
incident from the right: # = * L + * R . 

We now express iJ c i_; on in terms of the fermionic and 
bosonic creation and annihilation operators |44|: 



el— ion — 



= y"dr^V/ cl -ion(r-R 4 ) 



(28) 
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FIG. 4: Molecular junction undergoing local heating as a re- 
sult of transport. The electron gas has a temperature T e , and 
the phonons a temperature T„ . 
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(29) 



2u 



30 



X^4r^<4 2 (^+^), (30) 



since we have (Q;) a = £ i/3evib . A i<*,30 (&i/J+&J/j)- 
The clectron-phonon coupling constant is given by 



jia,<px 

El ,i?2 



Jdrj dK,|*^(r,K||) 
x^^ ps (r,R?)^ 2 (r,Kn), (31) 



where V ps (r, R?) denotes the pseudopotential due to the 
« th atomic core. It is interesting to note that, unlike the 
equilibrium case, in a current-carrying wire the electron- 
phonon coupling constant depends on two types of sta- 
tionary states: left- and right-moving states. 

We now assume that after continous exchange of en- 
ergy, both the electronic and phonon subsystems have 
reached a steady-state temperature and therefore assign 
a temperature T e to the electron gas, and a temperature 
T u to the phonons (see schematic in figure 0J, so that 
the statistics of the (jP) th phonon in the junction can be 
described by the Bose-Einstein distribution: 



1 



9iP 



where T u is the local temperature of the junction. 



(32) 



We can now use the Fermi's Golden Rule to infer rates 
of phonon emission and absorption, and thus explicitly 
evaluate equation Accordingly, the com- 

posite electron-phonon system's transition rate from an 
initial state \i) to a final state |/) is 



Rjf, = 2n\(i\H^_ ioii \f)\ 2 5(E f -Ei- uj jP ). 



(33) 



Ei and Ef denote the energy of the electron in the initial 
and final state, respectively. We are only considering the 
transition rate for an electron interacting with the (j/3) th 
mode, and so we have defined Hi? ■ to be the (ji3) th 

} cl — ion w ' / 

term in the H e i^ lon Hamiltonian (see equation H30(l '). 

Let us consider the case where an electron incident 
from the right emits a phonon. We can evaluate the 
statistical average of the square of the matrix element of 
the Hamiltonian in equation (|33|l to obtain 



(i\H 



E 



cl— ion 



/) ) 



(34) 



2uj 



1 A t2Q,LR 

Ji ia,j/3 J E-uj iB 



J/3 



:V /(l+^)/f(l-/^) 



Note that the 1 that is added to g^p corresponds to spon- 
taneous emission. 

The power emitted by right-incident electrons is sim- 
ply 

R,2 



w 



30 



■(Ei,E f ) = (E f -Ei)R jf3 (35) 

= ^\miLoJf)\ 2 (Ef-Ei) 

xSiEf-E.-cujp). (36) 

The total power emitted by right-incident electrons to the 
(j(3) th mode is the sum over all initial and final states: 



^ 2 = 2EE^ 2 (^) 



Ei E f 

2^£2n\(i\H. 

Ei E f 

xS{E f -E t 



l — i. 



Uj0) 



2 JdEiD%JdE f D\2-K\{i\H i —ion 
xiEf-EiWEf-Ek-Wjf,) 



f)\ 



(37) 



(38) 



(39) 



[i\H el _ ion \f)\ w jV j,(40) 



where we have multiplied by a factor of 2 due to spin 
degeneracy. In taking the continuum limit, we have in- 
troduced Z?£^ , which is the partial density of states of 
electrons moving to the right (left) with energy E. 

Finally, the total power emitted by right-incident elec- 
trons to the (j(3) th mode is therefore 



{Wff) = 2w(l+g jf) ) JdE X A ~>; 

lot 



ia,LR 



(41) 
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FIG. 5: Absolute magnitude of the power due to electron- 
phonon interactions as a function of T u - The intersection 
of the two curves gives the steady-state temperature of the 
junction. 




Bias (mV) 

FIG. 6: Steady-state temperature as a function of applied bias 
for two different contact geometries: a molecular junction, 
and a single gold atom contact. A schematic of the contact 
geometry is shown as an inset for each case. For both contact 
geometries, the nanojunction is not in thermal contact with 
the bulk electrodes. Reprinted with permission from 44]. 

Similar considerations apply to all other processes that 
contribute to equation Q22JI. 

The local temperature of the junction T w is evaluated 
(at a fixed electronic temperature T e ) when the total 
power (1221) is zero. This is illustrated in Figure where 
the magnitude of the sum of the heating processes is plot- 
ted on the same graph as the magnitude of the sum of the 
cooling processes as a function of T w . The point where 
the two curves intersect is the steady-state temperature 
of the nanojunction. 

Let us now discuss two examples. We first neglect cool- 
ing due to dissipation into the electrodes. Figure shows 




Bias (V) 

FIG. 7: Steady-state temperature as a function of applied bias 
for the same geometries as in Figure |H| where the dashed line 
corresponds to the molecular junction, and the solid line cor- 
responds to the gold point contact. Dissipation into the bulk 
electrodes is taken into account. Reprinted with permission 
from |44|. 

the steady-state temperature as a function of applied bias 
for two different nano junctions, a gold point contact and 
a molecular junction. We note that local heating occurs 
when a certain threshold bias is reached. This is because 
a given vibrational mode cannot be excited unless the 
incident electron is energetic enough to supply the requi- 
site amount of energy; that is, V"g™ = min{uj ia }. The 
onset bias for heating is in good agreement with exper- 
imental observations in the gold point-contact case [49(. 
Once the threshold bias is reached, however, the local 
temperature of the junction rises quickly, since cooling 
processes cannot effectively compensate for heating pro- 
cesses. Note that, in this case, very large temperatures 
can be reached at very small biases. 

This is in contrast to Figure [7| which depicts the same 
process, but with an important difference: the nano- 
junction and contacts are now assumed to be thermally 
coupled with each other, i.e. energy can be dissipated 
into the bulk electrodes. In order to estimate this en- 
ergy transfer we use the following expression for elastic 
phonon scattering between a bulk material A in contact 
with another bulk material B via a weak mechanical link 
(the nanojunction) |50| : 

I th =AirK 2 JdeeN A (e)N B (e)[g{T A ,e)-g(T B ,e)l (42) 

where g(T A rg\ ,e) is the Bose-Einstein distribution at a 
temperature ?a(b) and energy e and Na(t*) is the phonon 
spectral density of states of surface A(B) The weak 
mechanical link is modeled via a harmonic oscillator with 
stiffness K. 

Applying equation (|42|) to our case we obtain the re- 
sults of Figure [7] 44]. Notice that the bias scale in Fig- 
ure is much different than the one in Figure i.e. the 
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FIG. 8: Magnitude of the differential conductance through 
a molecular junction made of a benzene-dithiolate molecule. 
The dashed line gives the derivative of the differential con- 
ductance, while the insets illustrate the major long itudinal 
phonon modes Reprinted with permission from |§4']. Copy- 
right 2004 American Chemical Society. 



majority of heat generated in the junction is dissipated 
into the bulk electrodes. However, in real devices poor 
thermal contact between the junction and the electrodes 
can actually occur, as can localized phonon modes within 
the junction which have low coupling with the continuum 
of modes of the bulk electrodes |46( • Such instances can 
lead to very large local temperatures in the junction with 
consequent structural instabilities |44L l46l l51| . Such in- 
stabilities have been actually observed in metallic point 
contacts |52j and may be responsible for the low yield in 
fabricating molecular junctions 



INELASTIC CONDUCTANCE 

In addition to local heating, electron-phonon coupling 
can yield a lot of information on the internal structure of 
nanoscale junctions: discontinuities in the conductance 
can occur when the energy of incident electrons becomes 
large enough to excite different vibrational modes of the 
junction. 

In order to calculate the inelastic current, let us con- 
sider the case of an electron incident from the right. 
By treating the electron-phonon interaction as a per- 
turbation, the total wavefunctions of the (electron plus 
phonon) system can be written in terms of the states 
*B; n ict) = <8> \n ia ). The wavefunction of the sys- 
tem, including electron-phonon interactions is therefore 



(to first order): 



R., 
E- ' 



where 



\5tyg,ni a } is the leading-order change in the wavefunc- 
tion due to electron-phonon interaction. We can then 
calculate the leading correction to the total wavefunc- 
tion: 



lim > > 



dE'Dt 



— i> 



|*&n j7 J>l*&W'> 



e(E,n jl3 ) -e(E', n r >) 



Once again, D^ R ' is the partial density of states for 
electrons incident from the left (right). e(E, n ^) = 
E + (rijp + is the energy of state \^%; nj^.jsjal 

Using the identity lim e ^o jzj- — P(\) + i^8(z), where 
P(x) denotes the principle value of x, together with the 
matrix element that leads to equation l|34|) . we obtain 



IE 



(43) 



\^\- n jfl ) = ittJ2 J ^~ p A ^P De+^ y/mf^-fh+^UE^ \^e + ^e, n if > - 1) 



(44) 



The analagous expression for \8^>\\njp) is obtained by 
simply interchanging all R's with L's. 

As an axample let us assume that the electronic tem- 
perature T e is equal to zero. In that case, for an external 
bias V = E-ph — -^fr, only normal modes with energies 
ujjp < V can be excited. (Again, we are assuming the left 
electrode is positively biased.) Furthermore, if we assume 
negligible local heating, and that the few excited phonons 



decay on a short time scale, then we have gjp = 0. In 
this case, the inelastic contribution to the current is: 

pEpR p p 

SI = -i dE / dR / dK|, 

[{5*lyd z 8*l - d z (89%)*59%] . (45) 

In Figure we plot the inelastic conductance of a 
molecular junction. It is evident that only specific modes 
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with large longitudinal component (with respect to the 
direction of current flow) dominate the inelastic conduc- 
tance, while modes with large transverse component con- 
tribute negligibly. In addition, it is shown in Ref. p54| 
that inelastic current-voltage characteristics are quite 
sensitive to the structure of the contact between the 
molecule and the electrodes thus providing a powerful 
tool to extract the bonding geometry in molecular wires. 



CONCLUSIONS 

We have presented a review of several current-induced 
effects in nanoscale conductors and their description 
at the atomic level. These effects provide a wealth of 
information on the transport properties of atomic and 
molecular junctions beyond the value of the average 
current. In addition, their understanding is paramount 
to the possible application of nanoscale systems in 
electronics. 
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